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The collection efficiency in solar cells is treated by a new method in which 
all the effects of the solar spectrum and the absorption curve are contained 
in a single function readily obtained by numerical integration. The method 
is illustrated by a detailed study of the effects of surface recombination, body 
recombination and junction depth in silicon cells. The method is also gener- 
alized to include built-in electric fields, and calculations are given for sili- 
con. Sufficiently strong fields to improve the collection efficiency markedly 
can be produced in some compound semiconductors from a gradient in the 
energy gap. A discussion is given of the dependence of the collection effi- 
ciency on the absorption curve of the semiconductor. It is shown that silicon 
has a very favorable absorption curve in comparison with GaAs or InP. 
Finally, a treatment is given of the minority carrier collection in a two-junc- 
tion cell, and calculations are presented for silicon. It is concluded that this 
structure may be important for cells with high energy gaps and short life- 
times. 

I. INTRODUCTION 

The considerations reported in this paper have been stimulated by the 
current interest in the solar battery as a power supply for instruments 
and transmitters in satellite and space probe vehicles. A number of 
space vehicles t have contained solar batteries with peak outputs rang- 
ing from several watts to several hundred watts. It has been demon- 
strated that solar battery power supplies are technically feasible, not 
only in space vehicles but also in terrestrial telephone systems. 2 For the 
latter type of application, however, the solar battery has been found to 
be not competitive economically with several other available power sup- 
plies. 2 However, it is highly advantageous for space vehicles because of 
its advanced development and commercial availability, light weight, 
reliability and long life. The long life is due not only to the ruggedness 
and permanence of its structure and the absence of moving parts or 

t For a review of the use of solar batteries in space see Daniels. 1 
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chemically unstable components, but also to its external and independent 
source of energy, the sun. The belts of particle radiation 3 recently dis- 
covered surrounding the earth may materially reduce the life of solar 
batteries carried into space. ' We shall proceed with our discussion, 
however, on the assumption that the solar battery can survive the radia- 
tion. 

The solar battery is an array of hundreds or thousands of individual 
cells called solar cells. Each cell is a semiconductor slab, typically 1X2 
cm, containing a p-n junction within a very small distance, typically 
2 X 10 -4 cm, of the illuminated surface. The front and back surfaces are 
fitted with "ohmic" contacts for making electrical connection, and the 
front surface may be specially treated to reduce its reflectivity. 6 Com- 
mercially available cells are made of silicon, and have efficiencies of up 
to 14 per cent for converting the solar radiation incident upon them into 
electric power. 7 Other methods now known for converting solar radiation 
are far less efficient. f Thermoelectric converters, for example, on which 
considerable work has been done, can approach efficiencies of 1 per cent. 

The solar cell was invented by Chapin, Fuller, and Pearson, 10 who 
briefly described its fabrication, its principles of operation, and the limi- 
tations on its efficiency. Due to the work of the inventors and later au- 
thors, notably Prince, 11 Pfann and Van Roosbroeck, 12 Cumerow, 13 Ritt- 
ner 14 and Loferski, 16 the solar battery is well understood in terms of 
concepts familiar in electrical circuits and semiconductor physics. 

High efficiency in a solar battery would be desirable in any application, 
but especially so in space vehicles. The value of the vehicle, launched at 
considerable expense, depends in large measure on the instruments and 
transmitters it carries and the power available for this equipment. At 
the same time, the more equipment that is carried, the less space and 
weight can be allowed for the power supply. Therefore, in the economy 
of space vehicles the solar battery should have the maximum possible 
efficiency irrespective of the costliness of the improvements. 

The operation of the solar cell and the losses of efficiency can best be 
described in five steps: 

(a) Radiation is incident upon the surface and some is reflected with- 
out entering the cell. This reflection can be a very important loss, since 
the reflectivity 16 of clean silicon is about 30 per cent in the wavelength 
range of interest (0.4 to 1 n), and other materials that might be used 
also have high reflectivities. In practice it is found that the processing 
in the manufacture of silicon cells leaves the surface with quite a low 
reflectivity. 6 A recent study 16 has shown that the best antireflection 

t For a review of the utilization of solar energy, see Ref. 8, especially Vol. V. 
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treatments can increase the short circuit current of silicon solar cells by 
20 to 25 per cent compared with clean surface cells. It is reasonable to 
expect that the reflection loss from any material that might be used in 
the solar cell could be minimized in the same way. 

(b) The light enters the cell and some is absorbed by the intrinsic 
absorption process in which a hole-electron pair is created and a photon 
is destroyed.! The light absorbed in this way is the useful light in the 
solar cell. Light of wavelength longer than the intrinsic absorption edge 
cannot produce hole-electron pairs and is wasted in the solar cell. There 
is a further waste of energy when hole-electron pairs are produced by 
photons with more than the minimum required energy, since the excess 
energy is transferred in a very short time to the semiconductor lattice 
in the form of heat. These losses may be ascribed to the spectrum of 
solar radiation, since they would not occur if the radiation were mono- 
chromatic at the wavelength of the absorption edge. For silicon the 
losses due to the solar spectrum are about 53 per cent of the energy which 
enters the cell. 15 By choosing a semiconductor with a somewhat higher 
energy gapf this loss can be considerably reduced. ' ' The energy gap 
is therefore very important in considering materials for the solar cell.§ 

(c) Some of the minority carriers produced by the light flow by dif- 
fusion to the p-n junction. These are the carriers which contribute to the 
output current of the cell. Other carriers diffuse away from the junction 
and recombine at the surface or deep inside the cell. The percentage of 
minority carriers which contribute to the current is called the collection 
efficiency. In a typical commercial silicon solar cell of 9 per cent over-all 
efficiency the collection efficiency is about 60 per cent. 18 Minority carrier 
recombination is, therefore, a serious loss of efficiency in the solar battery 
at the present time. From the spectral response 6 ' 10,18 of silicon solar cells 
we can infer that the surface recombination velocity is very high, proba- 
bly greater than 10 5 cm/second. From the analysis to be presented we 
can also infer that the body lifetime is about 10 microseconds, the dif- 
fusion length|| about 10~ 2 cm. In high-purity silicon the lifetime can be 
several milliseconds 20 ' 21 and with certain surface treatments the reeom- 



t For a review of the optical properties of semiconductors, see Hrostowski. 17 

\ For the purposes of this paper the energy gap and the intrinsic absorption 
edge are essentially the same thing, except that the former is expressed as an 
energy while the latter is the equivalent photon wavelength. 

§ We shall not consider in this paper composite cells such as those suggested 
by E. Jackson (Ref. 8, Vol. V), in which one attempts to reduce spectrum losses 
by stacking several thin cells of different energy gaps. Cells of this type appear 
to be somewhat impractical from the standpoints of mechanical construction and 
of providing a suitably matched electrical load. 

fl Recently a measured value of ~5 X 10 -3 cm for the diffusion length has been 
reported by Vavilov, Smirnov and Patskevitch. 19 
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bination velocity 22 can be as low as 40 cm/second. The recombination, 
therefore, is due almost entirely to the degradation 2 of surface and body 
lifetime that occurs in the manufacturing process. In compound semi- 
conductors the recombination losses are likely to be considerably greater 
than in silicon. In gallium arsenide, for example, the lifetime seems at 
present to be of the order of a millimicrosecond. In general, it appears 
difficult to obtain lifetimes greater than 0.01 or 0.1 microsecond in the 
compound semiconductors, f 

(d) The diffusion maintains an excess concentration of minority car- 
riers on both sides of the junction. The voltage developed by the solar 
cell is due to these excess concentrations of minority carriers. This volt- 
age, however, is considerably less than the energy (in units of electron- 
volts) of a hole-electron pair in the semiconductor. The latter, for our 
purposes, may be taken to be the energy gap, which in silicon is 1.2 
volts. 25 The voltage of a silicon solar cell in full sunlight under maximum 
power conditions is about 0.4 volts. 10 ' 18 Therefore, the cell is able to con- 
vert only a portion of the energy stored as hole-electron pairs into elec- 
trical work. The loss may be referred to as the junction loss. The junction 
loss should vanish and the voltage should approach the energy gap when 
the minority carrier density approaches the majority carrier density, a 
limit corresponding to infinite light intensity. In the other limit of zero 
light intensity, the junction loss causes the efficiency of the solar cell to 
approach zero. According to the equivalent circuit point of view of Pfann 
and Van Roosbroeck, 12 the short-circuit current of the cell flows partly 
through the load and partly through the junction in the forward direc- 
tion. The voltage and the junction loss therefore depend upon the for- 
ward current-voltage characteristic of the junction. The theory of p-n 
junctions 26,27 predicts that the forward current should decrease exponen- 
tially with increasing energy gap. Therefore, insofar as actual junctions 
obey the ideal junction theory, the junction loss can be reduced by in- 
creasing the energy gap. A number of authors ' ' have considered the 
spectrum loss and junction loss as a function of energy gap. If all other 
losses are neglected, the maximum efficiency is obtained for an energy 
gap of about 1.6 volts. 15 It is now possible with mixed semiconductors 2 ,28 
to obtain nearly any desired energy gap from 0.7 volt (germanium) to 
2.4 volts (GaP), which completely covers the range of interest. It should 
be kept in mind that silicon p-n junctions show a large contribution to 
the current from thermal generation and recombination through traps 
in the junction region. 29 The short lifetimes in the compound semicon- 
ductors suggest that trap effects may be even more important in those 

t Ref . 24, p. 58. 
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materials. Therefore, it may not be possible to obtain large reductions 
in the junction loss by increasing the energy gap. 

(e) Finally there is the loss due to resistance of the very thin side of 
the junction next to the surface and of the contact to the surface. This 
resistance places a lower limit on the depth of the junction. Prince has 
considered the optimum depth for silicon solar cells taking into account 
the resistance loss and the collection loss. In practice, the internal re- 
sistance of silicon cells" is between one and two ohms and the junction 
depth is between 1 X 10 -4 and 2 X 10 -4 cm. The elimination of resist- 
ance loss would increase the efficiency of a 9 per cent cell to about 11 or 
12 per cent. 18 Unlike the other losses considered here, the resistance loss 
is not characteristic of the material used in the solar cell, and the methods 
for reducing it will be similar for any material. 

It has been mentioned that the solar spectrum loss and the junction 
loss have been considered at some length. It is also well understood what 
to do about reflection and resistance. But the collection loss has not been 
thoroughly treated in the literature. Several calculations have been 
made using the approach of Pfann and Van Roosbroeck. 12 This may be 
called the monochromatic method. Light of a certain wavelength, and 
therefore having a certain absorption coefficient in the material, enters 
the solar cell. Solutions are obtained in terms of elementary functions 
for the minority carrier density, taking into account diffusion, surface 
and body recombination, the generation of carriers by the light and the 
boundary conditions at the junction and at infinite depth. The solution 
gives a certain diffusion current into the junction which is the short- 
circuit current of the cell. This current must then be averaged over the 
wavelengths in the solar spectrum to obtain the collection efficiency. The 
averaging requires a tedious numerical integration, because the mono- 
chromatic current is a complicated function of the absorption coefficient, 
which is in turn a rapidly varying (measured) function of wavelength. 
The method is sound and can account in a satisfactory way for the col- 
lection efficiencies observed in solar cells. 6 The method is not well suited, 
however, for a systematic discussion of collection efficiency and no such 
discussion has been given. Several authors ' ' have even argued that 
the collection efficiency can be considered unity in fundamental con- 
siderations on the solar battery, since lifetimes can be expected to in- 
crease as technology improves. A review of the history of the photovoltaic 
effect and its utilization has recently been given by Rappaport. 30 

In this paper we consider the collection efficiency more systematically 
by a more powerful method than the monochromatic method. This 
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method is based upon obtaining a function, called the photodensity func- 
tion, which includes all the effects of the solar spectrum and the wave- 
length dependent absorption coefficient. This function is obtained by a 
relatively easy numerical integration over the solar spectrum. In terms 
of the photodensity function the solution can be obtained almost im- 
mediately to any solar cell collection problem in the approximation in 
which a single diffusion length describes the minority carriers. Illustra- 
tive calculations are presented for silicon solar cells showing how the 
collection efficiency depends on surface recombination velocity, junc- 
tion depth and diffusion length. The method is readily generalized to in- 
clude cases in which different diffusion lengths must be used for electrons 
and holes. Another generalization is presented which takes into account 
a "built-in" electric field in the region between the junction and the 
surface. It is shown that by the use of mixed semiconductors it should be 
possible to obtain sufficiently large built-in fields to increase the collec- 
tion efficiency significantly. In ordinary silicon cells, however, one would 
not expect the field to be large enough to have much effect. The im- 
portant question of which semiconductors should be best for solar bat- 
tery applications, already much discussed 11,13 ' 14, 15,3 ° with respect to spec- 
trum and junction losses, is taken up again from the point of view of the 
collection efficiency. It is pointed out that the absorption coefficient as a 
function of wavelength is very important in determining the collection 
efficiency. Silicon has an absorption curve of favorable shape, which in 
part accounts for its present superiority over other materials of more 
favorable energy gap. It is possible that for room-temperature use silicon 
will remain the best material, although higher-gap materials will cer- 
tainly be needed for use at temperatures above 200°C. Finally, there is 
presented a discussion of cells containing two junctions to improve the 
collection efficiency. It is shown that in cells of comparatively low col- 
lection efficiency ( <50 per cent), considerable improvement can be ob- 
tained by the use of a second junction. This construction may prove im- 
portant in high gap cells which might otherwise have rather low 
collection efficiencies. 

II. FORMULATION OF THE PROBLEM 

If the concentration of minority carriers is small compared to the con- 
centration of majority carriers, the equation describing the production, 
diffusion, and recombination of minority carriers is 

D <fn n + f° NMaMe -<" dx = o. (1) 

ax- t Jo 
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In (1), n(x) is the excess minority carrier concentration over the equi- 
librium concentration. In a typical case the maximum value of n(x) is 
of the order 10 12 cm -3 , so that the validity of (1) is assured. The diffu- 
sion constant D and lifetime t will be assumed to apply to all minority 
carriers whether holes on the n-side or electrons on the p-side of the 
junction. This assumption should lead to no serious error in silicon, but 
might have to be modified for some of the III-V semiconductors because 
of the relatively low mobility of the holes, f The integral term in (1) 
represents the production of minority carriers by light with a photon 
distribution N(\) in a material with absorption coefficient a(\) and 
intrinsic absorption edge X G . The total flux of photons capable of pro- 
ducing hole-electron pairs is 

AT = ( N(\) d\. (2) 



We shall neglect reflection completely and identify N(\) with the solar 
photon spectrum with respect to wavelength. J The total effective photon 
flux IsN = 3.3 X 10 17 cm -2 sec -1 for silicon.! The boundary conditions 
to be imposed on n(x) are 

n(«o)-0, n'(0) = («/2>)n(0) f (3) 

where s is the recombination velocity of the surface. A solution of (1) 
and (3) over the whole domain ^ x < °° represents the minority 
carrier density in a homogeneous illuminated semiconductor. 

If we now locate a junction at depth x = a, additional boundary con- 
ditions must be satisfied at the junction. In general, these will relate the 
minority carrier densities on each side of the junction to the operating 
voltage of the cell. The simplest case is the short-circuit condition, in 
which the voltage is zero.|| For this case the excess carrier densities must 
vanish on each side of the junction." 6,2 ' The boundary conditions for the 
short-circuit condition are therefore 

n(a) = 0. (4) 



tRef. 24, p. 12. 

j Details on the solar spectrum are given by Kef. 31. See also Ref. 15. 

§ This number, obtained from Ref. 31, is in substantial agreement with the 
plot of J* N(\)<1\ in Ref. 15. We are considering the solar radiation in space just 
outside the earth's atmosphere. 

|| We neglect voltage drops due to internal resistance. This causes no loss of 
generality, since an arbitrary resistance can be included in the equivalent circuit. 
In the presence of resistance, (4) corresponds to a small forward bias externally 
applied to the cell. The statement following (4) remains valid, and one identi6es 
this maximum diffusion current, not the short-circuit current, with the current 
generator. 
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This corresponds to the maximum current which can be drawn from the 
cell; smaller currents correspond to nonvanishing values of n(a) on both 
sides of the junction. However, we may regard any solution of (1) cor- 
responding to arbitrary operating conditions as the superposition of the 
short-circuit solution and an appropriate solution of the homogeneous 
equation 

D^- n = 0. (5) 

ax- t 

According to p-n junction theory, 26,27 these solutions are just those asso- 
ciated with forward currents in an unilluminated p-n junction. The super- 
position of solutions for the minority carrier density is therefore equiva- 
lent to superposing the short-circuit current with an appropriate for- 
ward current in the junction. If I a is the short-circuit current and Ij the 
forward current, the current in the external circuit is I g — I f , in ac- 
cordance with the equivalent circuit of Pfann and Van Roosbroeck. 12 

The collection problem, therefore, is to calculate the short-circuit solu- 
tion from (1) and (4) and from this the diffusion current density at the 
junction, taking into account the contributions from each side. The col- 
lection efficiency then is given by 

Q, = DN~\ | n'(a-) \ + | n'(o+) I ), (6) 

with the subscript "1" indicating that the expression applies to cells 
with a single junction. In a many junction cell Q n would be a sum of 
terms of the form (6). 

The equivalent circuit for a solar cell with load therefore consists of a 
current generator I g connected to a junction and a load in parallel, with 
the currents 7/ flowing in the junction (forward) and I f — I a in the 
load. Pfann and Van Roosbroeck 12 obtained the following condition for 
maximum power delivered to the load: 

G = z\nz + z - \, (7) 

where 



cV/kT 

z = e 



(8) 



V is the voltage, and G is a dimensionless reduced current proportional 
to I g . It is assumed that // obeys an ideal junction characteristic 26, 

I, = eAJ (e eV,kT - 1), (9) 
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where A is the junction area and J is a characteristic particle current 
density. The reduced current G is then defined by the relation 

G = I /cAJ = (N/Jo)Q. (10) 

The over-all efficiency can be written 

kT n (In s) 2 ni s 

e= F Q2 -^' (H) 

where z is the solution of (7), and 

W = AT 1 f N(\)(hc/\)d\. (12) 

Jo 

Clearly W is the radiant energy asborbed on the average to produce one 
hole-electron pair. For large G, (7) becomes 

Gttz\nz G»l, (13) 

and (11) becomes 

ett(eV/W)Q. (14) 

Here eV is the work done on the load and W/Qi is the energy absorbed 
per carrier flowing in the load. The approximation of large G is fully 
justified, since in full sunlight G will be of the order 10 7 . The efficiency 
for large G can also be written 

£l « (kT/W)Q l In G, (15) 

which exhibits the logarithmic dependence of the efficiency on the light 
intensity. It will be noted that the over-all efficiency d depends upon 
the collection efficiency Qi through the factor Q x and also through G. 

III. THE PHOTODENSITY METHOD 

The collection problem consists in solving (1) subject to (3) and (4). 
We begin by observing that the function 

r N(\)T T ^ m e-°*d\, (16) 

Jo 1 — a-L~ 

where L = (Dt) 5 is the diffusion length, is a particular solution of (1). 
It satisfies the boundary condition at infinity but none of the other 
boundary conditions. If ambiguity arises due to the pole ah = 1 oc- 
curring in the range of integration, the integral can be taken as a prin- 
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cipal value. The pole is removed if we combine (16) with an appropriate 
solution of the homogeneous equation (5), and write 

n{x) p = f Xg iV(X)r 1 a(X) „.„ (e~ ax - <T WL ) d\. (17) 

Jo 1 — orL- 

This function is a solution of (1) and satisfies 

w(oo) = n(0) = 0, (18) 

which is the form taken by (3) when the recombination velocity s is 
very large. Therefore, n(x) p represents the density of excess minority 
carriers in a semiconductor without junctions illuminated on a surface 
with fast recombination. The derivative of n(x) p at the surface is 

n'(0) p = f° N(\) -^7 dk. (19) 

Jo 1 + aL 

We call n(x) p the pholodensity . It can be readily evaluated by numerical 
integration (see Appendix), since the integrand is a simple expression 
with no poles. 

It is convenient now to introduce dimensionless quantities 



(20) 



f = x/L, (8 = aL, 7 = a/L, 
v(\) = N(\)/N, co = D/sL, 
and define the function 

<p(z) = (e z - l)/z. (21) 

Then the photodensity can be written 

n(x) p = (Nt/L)FU,L), 
where F(£,L) is the photodensity function 

F(t,L) = f e - r I ° v(X) 3-^^(1 - 0)] d\. (22) 

This function is shown in Fig. 1 as a continuous function of f for three 
values of L (10~ , 10 -3 , 10 -2 cm) based on the solar spectrum and the 
absorption curve for silicon. 3 "'! The calculation of F (f,L) is discussed 
in the Appendix. For small f the photodensity function has the ex- 
pansion 

HM = r f G "(A) /-^ dX - *** f "to***) d\ + ---. (23) 
Jo l + p Jq 



t Values of a > 10 6 cm • required as discussed in the Appendix for the cal- 
culation of Fig. 5 are from Pfestorf . 33 
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Fig. 1 — The photodensity function (22) for silicon as a function of f for L = 
l(r 4 , lO -3 and KT 2 cm. 

Once F({,L) 1ms been obtained, the solution n(x) satisfying (1), (3) 
and (4) is readily constructed from F({,L) and the functions e ±r , which 
satisfy the homogeneous equation (5). The solution may be written 



n(x) = (Nr/L)[Ae { + Be r + F(f,L)], 

where for ^ x £ a 

. 1 F(X)(1 + co) + co<r 7 F'(0) 



sinh 7 + co cosh 7 
= 1 F(y)(l - &,) + coe 7 F'(0) 



(24) 



(25) 



and for a ^ .r < °o 



sinh 7 + co cosh 7 

.4 =0, 

B = -e->F(y), 



(26) 



where F' = dF/d{ and for brevity L has been omitted as an argument 
of F. The solution in the region a ^ x < =0 is, of course, independent of 
co. In the region ^ .< ^ a the solutions in the limiting cases co — > 
and co — ► x are 



u -.o L L - smn y 



]• 



(27) 



„( X ) _ V* \ m _ F ( 7 ) + **h-r) ^(0)1 . ( 28 ) 
U ^.M L I cosh 7 J 
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The boundary conditions satisfied by these solutions at x = are 
n(0) = and n'(0) = respectively. 

The collection efficiency for single junction cells can now be obtained 
from (6) in the following form: 

Q l = ^(-r)(l + 0)^ + 0^(0) _ (29) 

sinh 7 + w cosh y 

In the limits o> — > and o> — > w this reduces to 

Q x -> F( 7 )(l + coth 7 ), 

W "° (30) 

Qi -* ^(7)(l + tanhi) + F'{0) sech y, 

Ci)-»0O 

which refer to fast and slow surface recombination respectively. 

Much of the simplicity of (29) and (30) results from the assumption 
of "effective" values for D and t which apply to both sides of the junc- 
tion. The difficulty of measuring these parameters in an actual cell, es- 
pecially as a function of depth, justifies this simplification for most con- 
siderations. The most conspicuous case where this assumption may lead 
to serious error is in some of the compound semiconductors where the 
mobility and diffusion constant may be an order of magnitude less for 
holes than for electrons, f When different L obtain on the two sides of the 
junction (24), (25) and (26) are still formally valid but (29) is no longer 
valid. The collection efficiency can then be written Qi = Q_ + Q + , 
where 

_ /''(TKcosh 7 + a> sinh 7) + coF'(O) „,, . 

Q- = T-r : -. t Ky), 

sinh 7 + o> cosh 7 (31) 

Q + = F{y) + F'{y) 

are the contributions from ^ x ^ a and a 25 x < °o respectively. If 
the same L is used for both sides of the junction Q- + Q+ reduces to 
(29). 

The fast surface recombination limit o> = is valid if w satisfies the 
two conditions o> « 1, o> <K 7. In sillicon solar cells the junction depth is 
considerably less than the diffusion length (7 « 1) so that the second 
condition implies the first. Therefore the criterion for high surface re- 
combination velocity is 

s » D/a. (32) 

The diffusion constant D can be obtained from measured mobilities n by 

t Ref. 24, p. 12. 
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means of the Einstein relation, 34 which in conventional laboratory units 



is 



D (in cm 2 /second) = 0.026 (~qqoj?) m (in cm 2 /volt-second), (33) 

with T the absolute temperature in °K. The mobility in highly doped 
silicon has been studied by Backenstoss, 35 who finds that for impurity 
concentrations greater than 10 19 cm -3 the mobilities of electrons and 
holes approach the limiting values of 80 and 40 cm 2 /volt-second respec- 
tively. Since a ~ 10 -4 cm, the criterion (32) gives s > 10 5 cm/second in 
silicon. Although surface treatments are known 21 which can reduce s in 
pure silicon to below 10 2 cm/second, Prince and Wolf 6 point out that 
these treatments tend to increase significantly the reflectivity of the 
surface. If this occurred the treatment might well reduce rather than 
increase the over-all efficiency of the cell. On the other hand, Malovet- 
skaya, et al. 16 have reported antireflection films that greatly reduce re- 
flection losses without increasing the surface recombination. The effect 
of surface recombination on Qi is shown in Fig. 2 for a silicon solar cell 
having the typical parameters 



D = 5 cm 2 /second, 
r = 2 X 10~ 5 second, 
L = 10 _2 cm, 
a = 10 _ cm. 



(34; 
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Fig. 2 — The collection efficiency as a function of surface recombination ve- 
locity for a silicon cell having the parameters (34). 
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For this cell Qi increases from 63 to 87 per cent as s decreases from 10 7 
to 10 3 cm/second. This example shows that considerable improvement in 
Qi would result from the elimination (s < 10 3 cm/second) of surface 
recombination. The presence of surface recombination in commercial 
silicon cells is indicated by the rapid fall-off in the spectral response 10,18 
at short wavelength X < 6 /z. Little is known about surface recombina- 
tion in the compound semiconductors. The fact that body lifetimes tend 
to be short 24 (t ~ 10~ 8 second) suggests that surface recombination is 
probably very fast (s ~ 10 cm/second). On the basis of these considera- 
tions, the most reliable guess is to set to = except where specific infor- 
mation to the contrary is available. 

The dependence of Qi upon the junction depth is shown in Fig. 3 for 
a silicon cell with D, t and L as given in (34). The curves shown for 
w = and w — > oo are envelopes for the entire family Qi(y,a)). Although 
the limiting curves have quite different forms they approach the same 
form for small and large depths: 



Ql ^ F'(0) = / \(X)/J(X)(1 -r-^r^X, 

•y-*0 Jo 



Q, -+ 2F(y). 

7 -»oo 



(35) 



The monotonic decrease of Qi(7,0) with increasing y is typical of silicon 
cells. In semiconductors of considerably higher energy gap the condition 
F'(0) + \F"{0) > may be satisfied, in which case Qi(y,0) will at first 
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Fig. 3 — The collection efficiency as a function of junction depth for a silicon 
cell with parameters I),t,L given by (34). The abscissa is y = a/L, where a is the 
junction depth and L the diffusion length. The curve Qi(t,°o) refers to zero and 
Qi(yfl) to infinite surface recombination velocity. 
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increase with 7 to reach a maximum and then fall to zero as 7 is increased 
further. A very broad maximum in Qi(7,«> ) is seen in Fig. 3 near 7 = 
0.2. The shape of the curve shows that relatively deep junctions (7 ~ 
0.5) could be used if surface recombination were absent, which would 
permit the internal resistance of the cell to be greatly reduced. On the 
other hand, deep junctions cannot be used if the surface recombination 
is important. This calculation also shows that in a cell with surface re- 
combination reducing a from 10~ 4 (7 = 0.01) to 5 X 10~ 8 cm (7 = 0.005) 
would improve the collection efficiency only slightly. Therefore, one may 
regard a = 10~ 4 cm as close to the optimum junction depth for silicon 
solar cells. 

The dependence of Qi on diffusion length is shown in Fig. 4 for a sili- 
con cell with a junction depth a = 10 -4 cm. The curves Qi(L,<x> ) and 
Qi(L,0) are envelopes for the family Qi(L,u) considered as continuous 
functions of L. These calculations show that L > 0.1 cm may be con- 
sidered essentially infinite. The limiting values as L — > =0 are 



Q X (L, 00) -* 1, 



\g 



1 - e 



(30) 



Qi(L,0) -> v(\) 

Jo <xa 



dX, 



which are indicated in the figure. It will be observed that Qi increases 
relatively little as L increases beyond 10 -2 cm. On the other hand con- 
siderable improvement results from increasing L from 10 -3 to 10 -2 cm. 
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Fig. 4 — The collection efficiency as a function of diffusion length for a silicon 
cell with the junction depth 10~ 4 cm. The curve Qi(L,=o) refers to zero and Q\(L,0) 
to infinite surface recombination velocity. 
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IV. THE EFFECT OF AN ELECTRIC FIELD 

In his review, Rappaport 30 has suggested that a "built-in" electric 
field may exist near the surface of solar cells causing minority carriers 
to drift toward the junction. If sufficiently strong this electric field would 
in effect eliminate surface recombination. One would expect that to re- 
duce surface recombination significantly the field would have to be of 
the order E ~ s/n ~ 10 3 volts/cm, assuming the typical values s~2X 
10 5 cm/second, y. ~ 200 cm 2 / volt-second. The built-in electric field could 
arise in two ways, from a gradient in the impurity concentration, and 
from a gradient in the energy gap. In the first case the field is given by 36 

E ~ (0.026/a) In (N /N a ) volts/cm, (37) 

where No ,N a are the carrier concentrations at x = and x — a respec- 
tively, and 0.026 = kT/e at room temperature. This is just the field 
required to cancel the diffusion current in equilibrium according to the 
Einstein relation (33). Actually (37) is an approximation giving an 
average effective field over the region ^ x ^ a. Consider as an example 
a cell made by the diffusion of boron from the vapor phase into w-type 
silicon. 6 ' 10 ' 11 The surface concentration of boron is about No ~ 18 cm" 
and N a ~ 10 17 cm -3 . If a ~ 10~ 4 cm, the built-in field is E ~ 500 volts/ 
cm, which is probably much less than s/n and therefore not large enough 
to cause much reduction in surface recombination. We may conclude 
that in commercial solar cells the built-in field can be neglected, as in 
the treatment of the last section and in the previous literature. On the 
other hand, the possibility remains that much larger built-in fields could 
be obtained, since JVo might be made to approach the solid solubility, 37 
which for boron exceeds 10 20 cm -3 . The built-in field would then be 
about E ~ 1800 volts/cm, which might cause a significant improvement 
in collection efficiency. 

It appears that really large fields E ^> s/n might be obtained from 
gradients in the energy gap. Such a gradient could be obtained in GaAs 
by diffusing in phosphorous from the surface. Mixed crystals 38 f GaAs- 
GaP can exist in all proportions, and the energy gap varies as a function 
of composition from 1.4 (GaAs) to 2.4 (GaP). If the impurity concen- 
tration is approximately constant, the band edge corresponding to the 
majority carriers must also be approximately constant. The gradient of 
the energy gap therefore appears almost entirely in the band edge for 
the minority carriers, and is equivalent to a built-in electric field acting 



t See also Ref. 24, p. 52. 
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on the minority carriers. It seems quite feasible that the field may be of 
order E ~ 10 4 volts/cm in GaAs-P. In addition to producing the field, 
the gradient in energy gap causes the absorption coefficient to become a 
function of position a(\,x). To make the analysis more tractable we may 
regard a(X) as an "effective" absorption coefficient independent of posi- 
tion. 

To take into account a constant electric field it is only necessary to 
add to the left side of (1) the term -nE(dn/dx), giving 

D^-"- liE ^ + f° N(\)a(\)e- az dX = 0. (38) 

dx 2 t dx Jo 

This is conveniently written in the reduced form 

^ _ y _ 28 | + f X ° tfe-K dX = 0, (39) 

af- as ■'o 

where 

E 



£ = 



V = 



2D/fxL ' 

n 



(40) 



Nt/L' 

and p,v£ are defined in (20). The boundary conditions on ?/(f) are 

y'(0) - Sy(0), y{y) =0, y(co) = 0, (41) 

where 

S = 28 + w _1 (42) 

depends on the electric field as well as the surface recombination ve- 
locity. If the last term in (39) is dropped the resulting homogeneous 
equation has the solution e~ pf , e^, where 

p = (1 + s 2 )*- Sgl, 

a = (1 + 8 2 )5 + £ £ 1, (43) 

pa = 1. 
The general solution can be written 

y(r) = Be-"" + Ae ai + F(f,L,8), (44) 
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where 

f X ° 8 e~ & * — e~ pX 

«'*»>- jr * il*\-U A (45) 

For the region ^ f ^ 7, 

F'(0) + (S - ff )F(7)e- 1 ' 



B = 



A = - 



p + o-e-^+'H + 5(1 - e -<"+*>?) " 



(46) 



As 8 approaches zero, F(£,L,8,) approaches F({,L) defined by (22) and 
A,B of (46) go over into A, B of (25). It will be assumed that the elec- 
tric field exists only in the region ^ f ^ 7. Thus the solution in the 
region £ ^ y is 

2/(f) = F(f,L) - F{y,L)e^. (47) 

Since // vanishes at the junction, the current is just the diffusion current 
and the collection efficiency is Qi = Q- + Q + , where 

Q + = y'(y + ), 

48) 
Q- = -Y(t-). 

From (47), 

Q + = F(y + ) + F'(y + ) (49) 

as in (31). From (44) and (46), 

_ F(y)[l - e- (p+,r)y + S(tr + p<r (p+ff)7 )] + F'(0)<r"(p + a) 



Q- 



p + < rc-c+ ff) > + £(1 - e -<^^) (5Q) 

- F'(y-), 



where F'(f) = (d/dS)F($,L,Z). In the limit 8 -> this goes over into 
Q_ given in (31). If 28 » 8, 28 » 1/8 and 28 » y, (50) becomes 

Qi -» [°v(l- e~ 0y ) d\-±-\° vSe-f dX. (51) 

Since the first term of (51) is the probability of absorbing a photon in 
the layer ^ f ^ 7, the limit 8 — > 00 corresponds to complete collec- 
tion of minority carriers produced in this layer. 
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A plot of Q- as a function of 8 is shown in Fig. 5 for a silicon cell 
with the parameters 

1) = 5 cm 2 /second, 

M = 190 cm 2 / volt-second, 



r = 



L = 



a = 



= 2 X 10 5 second, 



10 2 cm, 
10~ cm, 



(52) 



s = 2 X 10 5 cm/second, 

2D/nL = 5.3 volts/cm, 

s/n = 1050 volts/cm. 

The electric field is therefore E = 5.38 volts/cm. The value of Q+ cor- 
responding to these parameters is 0.52, so that the collection efficiency 
Qi = Q- + Q+ varies from 0.67 at zero field to 0.87 at infinite field. It will 
be observed that these values correspond in Fig. 2 to s = 2 X 10 5 cm/ 
second and s — > respectively, as one would expect. The critical field 
s/n indicated in the figure is the field at which about half the surface 
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Fig. 5 — The collection efficiency as a function of built-in electric field for a 
silicon cell with parameters (52). The abscissa is 8 = EnL/2D and 2D/ M L = 5.3 
volts /cm. The critical field 8/p and the asymptotic limits of Q\ are indicated by 
dotted lines. 
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recombination has been eliminated. Although it would be difficult to 
obtain built-in fields much larger than this in silicon, the high-held por- 
tion of the curve shows qualitatively the results to be expected from a 
gradient in the energy gap in a solar cell made of a compound semicon- 
ductor such as GaAs-P. The asymptotic formula (51) holds only for 
much higher fields than S = 10 4 because of the extremely high values 33 
of a near X ~ 0.3 n. The asymptotic limits at low and high fields are 
indicated by dotted lines. For estimating the effects of a given field it is 
convenient to regard the field as effectively lowering the surface recom- 
bination velocity. For estimating the effective velocity s* one may use 
the recipe 

8* = s[l + (En/s) (I + sa/D)]~\ (53) 

The collection efficiency can then be estimated from (29) by using an 
effective surface parameter w* = D/Ls* with considerably less labor 
than by evaluating (50) exactly. 

V. THE EFFECT OF THE ABSORPTION CURVE 

The absorption curve is the curve obtained by plotting the absorption 
coefficient a(X) as a function of wavelength X. Different semiconductors 
have quite different absorption curves as well as different energy gaps. 
In the several studies that have been made on the effect of the energy 
gap no consideration has been given to the absorption curve. 8 ' 11 ' 13,14 ' 16,30 
All of these studies have been concerned with spectrum and junction 
losses rather than collection losses. From the standpoint of collection 
efficiency, however, the absorption curve is of essential importance if 
there is any considerable surface recombination in the cell. Although the 
elimination of surface recombination may be in prospect, it is still of 
interest to consider the effect of the absorption curve in a cell with es- 
sentially infinite surface recombination velocity. This provides another 
standpoint from which to discuss the advantages of different semicon- 
ductors. 

In Fig. 6 are shown the absorption curves for silicon, 32 GaAsf and 
InP, 39 in the range 10 2 ^ a < 10 cm -1 . It will be observed that the curve 
for silicon differs from the other two in its gradual drop off with increas- 
ing wavelength. For a < 10 cm -1 it falls off more steeply, although not 
as steeply as GaAs or InP. The reciprocal a(X) _1 of the absorption coeffi- 

t The solid curve for GaAs is based on unpublished transmission measurements 
by W. G. Spitzer. The dashed part of the curve is an extrapolation based on a 
single additional reflection measurement by R. J. Archer giving a **** (1 ± 0.3) X 
10 6 cm- 1 at 0.546 p.. 
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Fig. G — The absorption curves for silicon, GaAs and InP. 

cient is a measure of the penetration distance of light of wavelength X. 
There is some wavelength Xi , for which 

a(Ai)" 1 ~ a, (54) 

which is the shortest wavelength of light useful in the cell. Light of 
shorter wavelength that Xi is absorbed too close to the surface to be col- 
lected at the junction. Similarly, there is a wavelength X 2 , for which 

o-CXs)- 1 ~ L, (55) 

which is the longest wavelength of useful light. If L = 10~ 2 cm or larger 
X2 is nearly equal to the intrinsic absorption edge X corresponding to 
the energy gap. Light of longer wavelength than X 2 is absorbed too far 
from the junction to be collected, or else is not absorbed at all. Since 
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the effective light lies in the range Xi ^ X ^ X 2 an effective photon flux 
can be defined : 



N ott = N(\) d\. 



(56) 



The collection efficiency may be written 

Qx = Nm/N, (57) 

where the flux N is defined by (2). The limits Xi , X 2 , and hence Q, evi- 
dently depend in an essential way on the absorption curve. 

The first equation (30) gives Qi for the present case in terms of the 
photodensity function F(^,L) denned by (22). Therefore, the qualita- 
tive considerations given above enter into determining the value of 
F(y,L). According to (22), F(y,L) contains an integral over X of the 
product of two functions 0/(1 + /3) and <p[y(l — (8)]. The weighting 
function v(\) representing the solar spectrum is slowly varying and not 
relevant to the present discussion. Fig. 7 shows plots of the two functions 
for a silicon cell with the parameters (34). It is clear from this figure 
that \i and X 2 correspond to cutoffs in <p and /3/( 1 + /3) respectively. The 
integrand in (22) is evidently small except in the region between the 
vertical lines in the figure corresponding toao = 1 and ah = 1. 
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Fig. 7 — The functions /3/(l + /3) and <p entering into the integrand of (22) 
defining the photodensity function. Calculation refers to a silicon cell with a = 
10 -4 , L = 10 -2 cm. Dashed lines show cutoffs Xi , X2 corresponding to (54) and (55). 
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It is now clear that when X x lies very close to X 2 the collection efficiency 
will be relatively low. This is the case with very steep absorption curves 
like that of In P. On the other hand, a gradual absorption curve like that 
of silicon leads to relatively high collection efficiency. If Xi , X 2 are de- 
fined by a(Xi) = 10 4 and a(X 2 ) = 10 2 cm -1 , (57) predicts the collection 
efficiencies 0.70, 0.21 and 0.045 for silicon, GaAs and InP respectively. 
The value for GaAs is subject to considerable uncertainty because it is 
based on a little-known part of the absorption curve. In any case, all of 
these numbers are qualitative, and mainly of interest to show that if 
surface recombination is high silicon is much superior to GaAs or InP. 
It is important to notice that this should be true in spite of the higher 
energy gap of the other two materials, because the differences in effi- 
ciency predicted for these materials on the basis of their energy gaps 
are much smaller than the differences in collection efficiency obtained 
here. For example, Loferski 15 finds that silicon and GaAs have relative 
efficiencies of 0.21 and 0.24 respectively, whereas silicon may have three 
times the collection efficiency of GaAs and 15 times that of InP. 

Several investigators have studied GaAs photocells. Gremmelmaier 40 
obtained an over-all efficiency of about 4 per cent and a collection effi- 
diency of about 0.2, in agreement with the estimates given above. An 
over-all efficiency of 6.5 per cent for a cell of very small area (0.007 
cm 2 ) has been reported by Jenny, Loferski and Rappaport. 41 These au- 
thors give current and voltage information for a somewhat larger cell 
(0.059 cm 2 ) having an over-all efficiency of 3.2 per cent, from which it 
can be deduced that the collection efficiency was about 0.26. It may be 
significant that the efficiency of 6.5 per cent was apparently obtained 
using a very low level of illumination (0.0057 watt/cm 2 ). Also, very low 
collection efficiencies have been reported for InP cells by Rappaport, 30 
although he attributes this to poor ohmic contacts to the cell. According 
to the present argument GaAs and InP will always give low collection 
efficiencies unless surface recombination can be eliminated, either by 
reducing the surface recombination velocity to below 10 4 cm/second, or 
by means of a built-in electric field of order 10 4 volts/cm. 

It should be mentioned that the spectral response for a GaAs cell re- 
ported by Jenny et al. 41 is not at all in agreement with the assumption 
of fast surface recombination. These authors observe an almost constant 
quantum efficiency from 5 to 9 y,. This differs radically from an earlier 
measurement by Seraphin,f which showed a pronounced fall-off at short 
wavelengths characteristic of surface recombination. No explanation has 
been advanced for this result nor can any be offered here. It may mean 

t Ref. 24, p. 65. 
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that surface recombination can be effectively eliminated in GaAs cells 
by treatments not yet defined, understood or controlled. 

VI. THE TWO-JUNCTION CELL 

In the typical silicon solar cell having the parameters (52), about 20 
per cent of the minority carriers produced are lost to surface recombina- 
tion and 13 per cent to body recombination. Therefore large reductions in 
surface recombination will be more rewarding in terms of over-all effi- 
ciency improvement than will proportionate reductions in body recom- 
bination. The two preceding secitons have been concerned primarily 
with surface recombination. It is also worthwhile to consider body re- 
combination in some detail, particularly since high-temperature applica- 
tions may be in prospect which will require high energy gap material. 42 
As discussed in the introduction, body lifetimes tend to be short in these 
materials. From Fig. 4 one may estimate that Qi is likely to be no larger 
than 0.2 in a cell with short lifetime. How this low collection efficiency 
may be improved by means of a second junction is discussed in this sec- 
tion. 

The collection efficiency for a two-junction cell will be denoted Q 2 . 
If the second junction is located at depth b below the surface, and 
i\ = b/L, the solution n(x) in the region a £ x £ b is 



, (58) 



L |_ sinh(i7 — y) sinh(ij — 7) 

which satisfies (1) and the boundary conditions n(a) = n(b) = 0. In 
the region x ^ b the solution is 

n(x) = Hi [F(f ) - e"- { F( v )}. (59) 

The solution in the region ^ x ^ a is the same as that obtained in 
Section III and given by (24) and (25). For the case of fast surface re- 
combination to = 0, this is 

satisfying n(0) = n(a) = 0. From (6), the total collection efficiency is 

Q 2 = F(y)[cothy + tanh $(»? - 7)] + F( v )[l 4- tanh \( v - 7)]. (61) 

As r] — > 7, Q 2 — > Qi for the case a> = given by (30). The effectiveness 
of the second junction may be measured by the quantiy 

«=«?!- Qi)/Qi , (62) 
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where Qi is the collection efficiency for a cell with a single junction at 
depth a. From (61), 

F(r,) 1 + tanh §(q - t) 1 - tanh $(* - t) 



5(tj,7,L) = 



F( 7 ) 



1 + coth 7 



1 + coth 7 



(63) 



Fig. 8 shows 8(r},y,L) plotted as a function of i\ for a silicon cell with 
a = 10 -4 cm and L = 10~ 4 , 10~ 3 and 10~ 2 cm. These calculations show 
that 5( tj) has a well-defined maximum which defines the optimum depth 
for the second junction for given a,L. 

Due to the loading effect of the second junction the relative improve- 
ment in over-all efficiency will be somewhat less than 5. This effect may 
be taken into account by use of the Pfann-Van Roosbroeck 12 efficiency 
expression (15). For the two-junction cell the reduced current may be 
written 



G 2 = Wi/Qx)Gx 



(64) 



where G\ = I„/eAJ = (N/J )Qi is the reduced current for the single- 
junction cell. Thus the efficiency expression for the two-junction cell is 



€ 2 « ■== Q2 In C7 2 
W 



(65) 



The relative improvement in over-all efficiency can be measured by 

5' = (e 2 - eOAi, (66) 
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Fig. 8 — The improvement of collection efficiency (defined by (62)) as a func- 
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where ei given by (15) is the efficiency of the single-junction cell with 
the junction at depth a. Finally, the over-all improvement 8' can be 
written as a function of the collection improvement 8: 



8' = 8 + (1 + 8) 



In j(l +_g) 
In (7 



(67) 



For practical purposes one may set In G = 18 in (67). The second term 
gives the effect of the loading due to the second junction. Two extreme 
cases may be noted: if 8 > 1, the second junction is a more effective 
collector than the first function and 5' > 8 ; if 8 < In J/ln G ~ 0.04, the 
second junction acts mainly as a load on the first junction and 5' < 0. 
A summary of calculations for silicon two-junction cells is presented in 
Table I. In each case it is assumed that the first junction is at a = 10 -4 
cm and the second junction at the optimum depth given in the third 
column. The collection improvement 5 and over-all improvement 8' are 
given in the fourth and fifth columns. In the second column is given the 
single junction collection efficiency in agreement with Q { (L,0) of Fig. 4. 
The last row with L = 10~ cm is the most typical of good silicon solar 
cells. For this case the improvement in collection efficiency is only 
5 = 0.092. On the other hand, for L = 10 -4 cm, 8 = 0.44 and 5' = 0.42. 
This case applies qualitatively to any short-lifetime cell, and shows that 
the two-junction structure may be useful for improving the efficiency of 
high energy gap cells. 



VII. SUMMARY 

In this paper the present status of solar cell theory has been reviewed, 
with emphasis on clearly defining the various mechanisms causing losses 
of efficiency. This review leads to the conclusion that the problem of the 
collection of minority carriers has not received attention in the literature 
commensurate with its importance. 

The collection problem is formulated and then solved by a new method 
in which all the effects of the solar spectrum and absorption curve are 
contained in a single function, the photodensity function. The method is 
convenient for most solar cell collection problems and especially so when 
a single diffusion length can be used. 

Table I 



L (in cm) 


Qi 


b (in cm) 


S (max) 


6' 


io-* 

10~ 3 
10~ 2 


0.186 
0.420 
0.634 


0.032 
0.013 
0.008 


0.443 
0.258 
0.092 


0.417 
0.226 
0.055 
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The familiar single-junction silicon cell is considered first and calcula- 
tions are presented to show how the collection efficiency varies with sur- 
face recombination velocity, junction depth and diffusion length. It is 
found that the elimination of surface recombination would not only im- 
prove collection efficiency but also permit the internal resistance to be 
greatly reduced. Little or no improvement in silicon cells is to be expected 
from making the junction depth less than 10~ 4 cm or the diffusion length 
longer than 10 -2 cm. 

The theory is extended to include a constant built-in electric field, and 
calculations are presented for silicon. It is concluded that commercial 
silicon cells do not have large enough built-in fields to affect the collec- 
tion efficiency significantly. Sufficiently large fields should be obtainable 
in some compound semiconductors from gradients in the energy gap. An 
approximate relation is given for an effective surface recombination 
velocity less than the true velocity, which takes into account the effect 
of the built-in field in reducing surface recombination. 

A discussion is given of the dependence of the collection efficiency on 
the absorption curve of the semiconductor. This discussion provides an- 
other basis beside the energy gap on which to compare different semi- 
conductors for solar battery use. It is shown that silicon has a very favor- 
able absorption curve in comparison with GaAs and InP. It is suggested 
that this accounts to a large extent for the continuing superiority of 
silicon over other materials with larger energy gaps. The theoretical 
superiority of GaAs and several other higher gap materials over silicon 
can only be realized if surface recombination can be drastically reduced ; 
specifically, it is necessary that the effective surface recombination ve- 
locity be reduced to 10 4 cm/second or lower. 

Finally, a two-junction cell is considered in connection with the prob- 
lem of reducing body recombination. Although body recombination is 
considerably less important than surface recombination in good silicon 
cells, it will be probably much more important in high energy gap cells. 
The illustrative calculations presented for silicon also have qualitative 
significance for other cells. It is shown that an improvement of about 42 
per cent in over-all efficiency may be expected from the two junction 
structure in material with a diffusion length of 10~ 4 cm. The improvement 
in a good silicon cell with a diffusion length of 10~ 2 cm would only be 6 
per cent. 
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Note. Since the completion of this work a review paper by Wolf 43 has 
appeared which discusses some of the topics taken up here. 



APPENDIX 

In this section some comments are given which may be helpful for the 
evaluation of the photodensity function (22). A very accurate and con- 
venient method is the Gauss quadrature, which approximates the in- 
tegral by a summation over the values of the integrand at certain specific 
(not equidistant) points. The limits 0,X o in (22) can be replaced by 
Ai , A2 so chosen that the contributions from the regions < A < \\ and 
X2 < X < X G are negligible. The approximation is to represent the integral 

F(f,L) = \° g(\)d\ 

Ja 



(08) 



10 

= Hh - Xj) T,KjgM- 



i=l 



The choice of a 10-point Gauss quadrature here is arbitrary, but has 
proved to be a satisfactory compromise between convenience and ac- 
curacy. For silicon one can choose Xi = 0.42 p, X 2 = 1.08 n, and (G8) 
becomes 



f g(\)d\ = £ (O.33# y )0(X;). 



(09) 



The wavelengths \j and corresponding absorption coefficients 32 a(Xy) and 
solar spectrum 31 weights v(\j) for this case are given in Table II, along 
with the quadrature weights 0.33ft/ . The integrand g(X) can be readily 
evaluated with a slide rule for each Xy . 

Table II 



Mm) 


Q.33R 


a (in cm ') 


Km" 1 ) 


0.429 


0.0220 


3.7 X 10 4 


1.22 


0.465 


0.0492 


2.0(4) 


1.50 


0.526 


0.0723 


9.0(3) 


1.57 


0.607 


0.0888 


4.3(3) 


1.62 


0.701 


0.0977 


2.2(3) 


1.51 


0.799 


0.0977 


1.03(3) 


1.36 


0.893 


0.0888 


4.5(2) 


1.26 


0.974 


0.0723 


1.56(2) 


1.12 


1.035 


0.0492 


42. 


1.05 


1.071 


0.0220 


17. 


1.02 
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The accuracy of (69) may be tested with an elementary integral hav- 
ing a qualitative behavior similar to g(\) in the range of integration. 
Such an integral is 

-1.08 

/ (1.08 - X)*(X - 0.42)* dX = 0.17100. 

JO. 42 

The quardature (69) with X, and 0.33A; from Table II gives 0.17111. 
The derivative F'($,L) is defined by the integral 

= r f f ° v{\) --^— Q {i - «vir(i - p)]\ rfx. 

Jo 1 "T P 

This integral may be accurately evaluated by (69), providing /3(Xi)f » I 
to insure that the truncation error at Xi is small. The photodensity func- 
tion for a constant electric field denned by (45) can be written 

F(f,L,8) = pf c- pf I ° v — £— *>[pf(l - ap)] dX, (71) 

Jo 1 ■+■ PP 

and its derivative is defined by 

F'(f,L,8) = pe-* 1° v — ^— {1 - ^[ P r(l - 00)]} d\. (72) 
Jo L -\r pp 

These expressions are in the form of (22) and (70) respectively, and 
can be evaluated by (69) in most cases. 

Some of the integrals occurring in the theory cannot be evaluated 
with sufficient accuracy by (69) because of large contributions from the 
short wavelengths X < Xi = 0.42 p. As an example of this difficulty con- 
sider 

I ° v(\)d\ = 1, (73) 

Jo 

which follows from the definition of v(\) in (20) and of N in (2). The 
quadrature (69) gives 0.905 for this integral, which shows that 

/ v f/X = 0.095. 

Jo 

This result can be used to correct (69) for the evaluation of 

F'(0,L) = f\^-d\ (74) 

Jo I "T P 
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occurring in (31). If /3(X X ) » 1 it is only necessary to add 0.095 to the 
right side of (69) ; thus for silicon 



F'(Q,L) = 0.095 + £ (0:MRj) Vj — ^_ 



(75) 



The evaluation of F'(0,L,&) occurring in (50) is more difficult and re- 
quires a numerical integration over the larger region < X < X G . To 
carry out the integration the absorption curve must be known for very 
high absorption 3 ' 1 a ~ 10 cm -1 . There may also be a significant correc- 
tion to be added to (09) in the evaluation of F(£,L,S) for very high fields 
S > 100. These corrections can be readily calculated by integrating 
from 0.22 to 0.42 n by the trapezoid rule. The necessary data for silicon 
is listed in Table III. For large fields (72) becomes 



F'(f,Z,,8)->;I f° V 0e-^ ( lX, 

f,-»=o -SO •'0 



(76) 



which can be evaluated by (69) if /3(Xi)f » 1. For silicon with a = 10 * 
cm the following values were found for the integrals appearing in (23), 
(36), (51) and (76): 



h 1 - 



/ 



— e 



aa 



padX = 3.7 X 10 cm" 



d\ = 0.754, 



/\(1 - c~ an ) dX = 0.351, 
J vae~ aa d\ = 1490 cm" 1 . 



(77) 



Table III 



X00 


r&r») 


a (cm-') 


0.22 


0.00477 


1.2 X 10 6 


0.2-1 


0.0148 


1.4 


0.2G 


0.0506 


1.5 


0.28 


0.0625 


1.6 


0.30 


0.178 


1.5 


0.32 


0.300 


1.3 


0.34 


0.448 


0.9 


0.36 


0.538 


0.4 


0.38 


0.654 


0.1 


0.40 


0.801 


0.06 


0.42 


1.14 


0.04 
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